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Abstract 

The time dependent neutrino oscillation signals due to the passage of a shock wave 
through the supernovae are analyzed for the case of three active neutrinos and also for the 
case that there are two additional sterile neutrinos. It is shown that, even without flavour 
identification and energy measurement, detailed information about the masses and mixing 
angles of the neutrinos may be obtained with a detector with excellent time resolution such 
as IceCube. Such a signal would also give important information about the nature of the 
shock wave within the supernovae. 



1 Introduction 

The possibility that, in addition to the three flavours of light neutrinos required by the Standard 
Model, there are light sterile neutrinos has been severely constrained by neutrino oscillation 
experiments. However there is still the outstanding problem of explaining the data from the 
Liquid Scintillator Neutrino Detector (LSND) |lj which, if confirmed, might require the existence 
of at least one extra light neutrino. The precision measurements at LEP have shown that only 
three neutrinos with standard weak interactions exist with mass less than and therefore 

additional neutrinos must be sterile 0. However the addition of a single sterile neutrino to 
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explain the LSND anomaly is not able to give a good fit to solar and KamLAND ^ ^ and 
atmospheric El data together with the data from the short-baseline experiments CHOOZ 
[71 Bugey, CCFR84, CDHS, KARMEN and NOMAD 0. While a 3+1 hierarchy in the four 
neutrino case is preferred over a 2+2 hierarchy, neither is acceptable [S| and one is driven to 
consider the case of two additional sterile neutrinos fU] • In this case a 3+2 mass hierarchy gives 
an acceptable fit, significantly better than that found for the 3+1 case. 

Pinning down the properties of additional sterile neutrinos is a daunting task. Even if the 
LSND result is not confirmed there will still be room for light sterile neutrinos mixing with the 
three active states that terrestrial experiments will not be able to exclude. Here we investigate 
the possibility of gaining significant new information on this sector from the neutrinos coming 
from supernovae. The reason this can give a significant signal is due to the possibility of resonant 
conversion of the various neutrino species (the Mikheyev-Smirnov-Wolfenstein (MSW) effect [TT] ) 
within the supernova. The MSW effect, being a resonant process, is much more sensitive to small 
mixing angles than non resonant oscillation phenomena. Feasibility studies of using neutrino 
signal from a galactic supernova to determine the neutrino mass hierarchy and constrain the 
mixing angle 613 have been done before in for the standard case of three active neutrinos 
and in jT^j for three active and one sterile neutrino. 

A further significant advantage of the supernovae signals for neutrino oscillation has been 
appreciated in recent years [14J. This follows from the fact that in the supernova explosion a 
shock wave propagates through the supernova. From the theoretical simulations of supernovae 
it is believed that when the core of the collapsing star reaches nuclear density, the collapse 
rebounds forming a strong outward shock. This is stalled and then regenerated by a neutrino 
driven wind. During this process it is believed that as well as a forward shock a reverse shock 
could also form At the shock front the density changes very rapidly so that the resonant 
transition may become non-adiabatic as the shock wave passes |14|, [T3]^. As a result one may have 
rapid changes in the active neutrino luminosity produced by the supernova, very characteristic 
signals which do not require knowledge of the overall luminosity. These signals are observable by 
detectors such as IceCube ^ capable of collecting high statistics with excellent time resolution. 
IceCube will provide a km^ neutrino detector with a time resolution of 10ns ^7]. Although 
designed to observe high energy neutrinos from astrophysical sources, it has been realized that 
the lower energy neutrinos from supernovae can be probed in IceCube through the detection of 
additional photons in the background halo, caused by supernovae neutrinos interacting with the 
ice [21] . IceCube cannot measure the energy spectra of the supernova neutrinos, but because of 
its excellent time resolution and high statistics, we find, for the reasons given above, that it can 
be an excellent probe of a sterile neutrino component to neutrino mass eigenstates. Moreover 
observation of such signals in effect provide information on supernovae seismology and can yield 
important information on the nature of the shock wave within the supernovae. 

For the case of oscillation between active (anti) neutrinos only, detection of such effects re- 
quires a difference in the initial properties of the neutrinos produced at the neutrinosphere. 
However there is considerable uncertainty in this. Early studies |22] suggested that, while the 
luminosities of the three antineutrino species should be quite close, the average energy of the z7e 

^While completing this paper we received a paper studying the effects of turbulence on the shock wave. 
These effects have not be included here but may be expected to broaden the structures presented below. 

^For detailed discussion of the effect of the shock wave on the signal in water Cerenkov detectors, we refer 
the readers to [14[I18[[T^ for effects of the forward shock only and |15[ I2U| for both forward and reverse shocks. 
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should be about half that of and Vt-. More recently studies [231 including additional scattering 
processes within the neutrinosphere have suggested that the difference in the average energies of 
the z/g and v^Ivt is quite small but that there is a significant difference in their luminosities. We 
will investigate the time dependent signals for active neutrino oscillation for both these cases. 
For the case of oscillation to sterile neutrinos the time dependent oscillation signals are not so 
sensitive to uncertainties in the initial antineutrino processes because sterile neutrinos produced 
by resonant oscillation will not be visible in detectors and so oscillation to a sterile neutrino 
is observable even with a detector having no capability to distinguish the flavour or energy of 
active neutrinos. 

In this paper we study the signals to be expected by IceCube both for the case of just 
three active neutrinos and for the case that there are two additional sterile neutrinos. To make 
the discussion tractable we concentrate on the range of mass and mixing angle parameters 
that provide an explanation of the LSND events while being consistent together with all other 
neutrino oscillation data. However our results indicate that such signals will also be significant 
in constraining the sterile neutrino component even if the LSND result is not confirmed. 

In Section 2 we briefly review resonant conversion within the supernova, the various possibil- 
ities for the neutrino spectra and list the parameters used in the subsequent analysis. In Section 
3 we discuss the properties of the neutrinos at the neutrinosphere and the nature of the shock- 
waves propagating through the supernova. Section 4 contains the results of our analysis. We 
first review the characteristics of the IceCube detector and then present its detection rates for a 
"standard" supernova for the case of three active neutrinos only and for the case of three active 
neutrinos plus two sterile neutrinos for the various possible neutrino mass hierarchy possibilities. 
Finally in Section 5 we present our conclusions. 

2 Resonant Conversion Within the Supernova 
2.1 Three Active Neutrinos 

We start with the case of three active (anti)neutrino species. In the basis of flavor eigenstates 
the evolution of (anti)neutrinos is described by the effective Hamiltonian 



where is a unitary matrix and is defined by |z/j) = Uic\vc). The vacuum mass eigenstates 
are Vi and the weak interaction eigenstates are v^. is the matrix of mass squares of the 
mass eigenstates in vacuo, and A is the matter induced mass matrix. For mixing between three 
active neutrinos 



n = ^{uM'\j^ ^A) , 




(2) 



A 



Dm(7(Ai,0,0) , 



(3) 



Ai = ±V2GFpNAYe X 2E , 



(4) 



3 




Figure 1: The evolution with density of the mass eigenvalues (squared) for the case of three 
active antineutrinos (a) for the normal hierarchy (Amgj^ > 0) and (b) for the inverse hierarchy 
(Am^i < 0). 

where i = 1 — 3 and a = e, /i or r, is the mass squared of the i*^ vacuum mass eigenstate, 
Gf is the Fermi coupling constant, p is the density of ambient matter, A^^ is the Avogadro's 
number, is the electron fraction and E is the energy of the neutrino. Ai is the matter potential 
induced by charged current interactions of Uf. (z/e) with electrons and the "+" ("— ") sign in Eq. 
(0} corresponds to neutrinos (antineutrinos). There is also an effective mass squared induced by 
neutral current interactions of the (anti)neutrinos with matter. However, this does not have a 
physical effect because the neutral current interaction is universal for all (anti)neutrino flavours, 
while oscillations depend on the difference of mass squared and not on the absolute value. The 
mixing matrix U is given in terms of mixing angles and CP-violating phases. If CP conservation 
is assumed the mixing matrix takes the form 

U= U fl^''^(^AB), (5) 

0<B<AA=2 

where R^^{6ab) are rotation matrices, representing a rotation of {6ab) in the AB plane. The 
mixing matrix for three flavours is therefore given as 

(C12C13 S12C13 Si3 \ 

-S12C23 - C12S23S13 C12C23 - S12S23S13 S23C13 , (6) 
~ Cl2C23'5i3 — C12S23 — S12C23S13 C23C13 J 

where Cij = cos 9ij and Sij = sin 9ij . 

The difference in mass squares of the mass eigenstates in matter can be calculated by diago- 
nalizing Eq. and is shown in Fig. ^as a function of the matter density for the antineutrinos^. 
Flavour oscillations predominantly occur at the resonance densities, where the effective mass 
difference between the two relevant mass eigenstates in matter becomes minimum. If the reso- 
nances are far apart then each resonance can be treated as an effective two neutrino problem, 

^Since the matter potential A is negative for the antineutrinos, we show the evolution of the mass squares of 
the mass eigenstates in Fig. QJas a function of — p, where p is the matter density. 



4 



independent of other resonances. For two flavor oscillations the resonance condition is given by 



Ai = Am^ cos 26 , (7) 

where Am^ = — ml is the mass squared difference and 6 is the mixing angle between the 
two states in vacuum. At the resonance a flip between the mass eigenstates is possible due to 
the changing mass density. The "flip probability" between the two mass eigenstates is given by 

exp(— 7 sin^ 6*) — exp(— 7) 



1 - exp(-7) 



7 = TT- 



Am^ 



E 



dlnAi 



dr 



(9) 



where r^m is the position of the maximum violation of adiabaticity (mva) |23j and is defined as 

Ai(r„,,,) = Am' (10) 

For a given set of values of Am' and 6, if the density gradient is small enough so that 7^1 
and 7sin'^ ^ 1, the flip probability reduces to Pj = exp(— 7 sin' 6') ^ and the resonance is 
called adiabatic. However, if for the given Am' and 6, the density gradient is large such that 
7 <C 1, then Pj ^ cos' 6. In that case for very small mixing angles Pj ~ 1 and the resonance is 
completely non-adiabatic. For all intermediate regions Pj ranges between [0 — 1]. 

For three active neutrinos with their vacuum masses given by the current experimental data 
and for a static supernova density profile, the resonance condition is satisfied at two distinct 
widely separated densities inside the supernova [221 ■ The resonance at the higher density is driven 
by Am']^ and ^13, while that at lower densities is driven mainly by Am'^^ and 612. Since the 
matter induced potential Ai is positive (negative) for neutrinos (antineutrinos), the resonance 
condition given by Eq. (|7j) is satisfied only if Am' is positive (negative). Since Amgi (for 
612 < 7i"/4) is known to be positive at a very high level of confidence from the solar neutrino 
data, the lower resonance is therefore satisfied only for the neutrino channel. However, Am'^ 
could be either positive (normal hierarchy) or negative (inverted hierarchy) and therefore the 
higher resonance can occur either in the neutrino or the antineutrino channel. It can be shown 
that for the most plausible supernova density profiles and the current allowed values of Am^i and 
612, the flip probability at the lower resonance is zero and the transition is completely adiabatic. 
However, depending on the value of ^13 and the neutrino mass hierarchy, the flip probability at 
the Am']^ driven higher resonance may have any value between [0 — 1] . We refer readers to |T2] 
for detailed discussion of the three generation oscillation probability in a supernova environment 
with static density profiles. 

Things get more involved when one considers the effect of shock waves on the supernova 
density profiles. As we will discuss in the following section, the effect of the shock is to cause 
very sharp jumps in the density gradient. This results in the same Am' producing multiple 
resonances which are relatively close together. If we assume that the phase effects can be 
neglected even in this case"^, then the individual resonances can be considered as independent 



■*This is a good approximation especially in our case since we will be working with the IceCube detector which 
does not have any energy sensitivity. This means that the detected supernova events will be averaged over energy 



and hence phase effects will be further washed out. 
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two generation resonances and the net flip probability Ph can be expressed in terms of the 
multiple flip probabilities Pi as |2ni UHj 

Ph 1~PhJ V l-P^J ' ^ ^ 

where n is the number of resonances occurring for the same Am^ due to the shock effect. 

2.2 Inclusion of Two Sterile Neutrinos 

For mixing between 3 active neutrinos and 2 sterile neutrinos the form of Eq. ^ remains the 
same but now we have 

f/ = n n i?^^ (12) 

B>AA=\ 

where R^^ is a 5 by 5 rotation matrix about the AB plane, 

= Diag{m\, ml, ml, ml, m\) (13) 

A = Diag{Ai,Q,{},A2,A2) (14) 
Ai = ±V2GFpNAYe X 2E (15) 

A2 = ±V2GfpNa{1 -Ye)xE (16) 

where, as in Eq. (@]), the "+" and "— " signs in Eq. (fT3j) corresponds to neutrinos and an- 
tineutrinos respectively. As before, i = 1 — 5 and a=e, p, r, si or S2, where Si and S2 are two 
sterile neutrinos. A2 = —A^c, where A^c is the matter- induced potential of the neutral current 
interaction of neutrinos/antineutrinos with matter. As discussed in the previous section, all the 
three active neutrinos/antineutrinos pick up an equal matter-induced potential A^c in matter 
and the matter induced mass matrix is given by 

A' = Diag{Ai + Anc,Anc,Anc, 0,0) , 

= A + Anc1. (17) 

Since the sterile neutrinos do not interact with matter, A^c cannot be completely factored 
out from the mass squared matrix as in the three neutrino case and hence can be recast as a 
negative matter potential for the sterile states. The effective mass squared differences between 
the neutrino/antineutrino eigenstates in matter therefore depend on both Ai and A2. 

The difference in mass squared can again be obtained by diagonalizing the mass squared 
matrix in matter and the results for the antineutrino channel are shown in Fig. |21for the different 
possible mass hierarchies of the flve (anti)neutrino system in vacuum. As discussed above, while 
we know that Am^i > 0, sign of Amg^ is presently unknown. The only experimental information 
we have on Amli and Am^i comes from the data of short baseline experiments including LSND 
and their sign are completely unknown as well. Therefore, in principle its possible to have six 
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Figure 2: The evolution with density of the mass eigenvalues (squared) for the case of three 
active anti-neutrinos and two sterile anti-neutrinos, for the mass hierarchies (a) N2-I-N3, (b) 
N2+I3, (c) H2+N3, (d) H2+I3, (e) I2+N3 and (f) 12+13. N3 and 13 correspond to the normal 
and inverse hierarchies in the active neutrino sector. N2 corresponds to the normal hierarchy in 
the sterile sector where Am^^ > and Amli > 0. H2 corresponds to the half hierarchy in the 
sterile sector where either Amg^^ > and Am|^ < 0, or Amg^^ < and Amli > 0. 12 corresponds 
to the inverse hierarchy in the sterile sector where Am|j^ < and Amli < 0. 
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types of mass hierarchy. For simphcity of notation, we will henceforth call them: 

(a) N2+N3: Am|^ > 0, Amj^ > and Am^^ > 0, 

(b) N2+I3: Am|i < 0, Amf^ > and Am^^ > 0, 

(c) H2+N3: Am^^ > 0, Amj^ > and Aml^ < 0, 

(d) H2+I3: Am|i < 0, Aml^ > and Aml^ < 0, 

(e) I2+N3: Am^^ > 0, Amj^ < and Am^^ < 0, 

(f) 12+13: Aml^ < 0, Am|i < and Am^^ < 0. 

As discussed in the previous section the passage of the shock wave causes each Am^ to have 
multiple resonances. The additional resonances due to the inclusion of sterile neutrinos can also 
be considered as independent, and the net flip probability is given by equation ITTl 

2.3 Experimental Bounds 

Here we very briefly review the current knowledge we have about the masses and mixing angles 
of the active and sterile neutrinos from all available neutrino oscillation data. 

2.3.1 Solar Data 

The combined analysis of world data on solar neutrinos [3, and the KamLAND reactor data j3] 
has established the so-called Large Mixing Angle (LMA) solution as the solution to the solar 
neutrino anomaly with the current best-fit parameters j2Zj AttIq = Am^i = 8.0 x 10~^ eV^ and 
sin^ 6*0 = sin^ 6*12 = 0.31 and the 3a allowed range given by j2Zl I2H] 

0.25 < sin^^ia < 0.39 (18) 

7.2 X 10-^eV^ < Aml-^ < 9.2 x lO'W (19) 

2.3.2 Atmospheric Data 

The zenith angle dependent event spectrum of the atmospheric neutrino data from Super- 
Kamiokande [5] and the data from the long baseline K2K experiment can be best explained in 
terms of almost pure — oscillations with best-fit parameters Aml^^ = |Am|]^| = 2.1 x 10~^ 
eV^ and sin^ 26'atm = sin^ 26*23 = 1 and the 3a allowed range given by [S] 

sin2 2^23 > 0.9 (20) 

1.3 X 10-^eV^ < Am|i < 4.2 x lO'^eV^ (21) 

2.3.3 CHOOZ Reactor Data 

The upper limit on the mixing angle ^13 is mainly determined by the reactor neutrino experiments 
CHOOZ and Palo Verde |7j. Data from these experiments when combined with the solar and 
atmospheric neutrino data gives at 3a the bound [2HI I2S] • 

sin^ Ois < 0.044 (22) 
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Parameter 


Best fit 


Best fit in sub-eV 




0.92 eV2 


0.46 eV^ 




22 eV^ 


0.89 eV^ 


Ue4 


0.121 


0.090 




0.204 


0.226 




0.036 


0.125 




0.224 


0.160 



Table 1: The two best fit solutions for the short basehne experiments taken from |10j . 
2.3.4 Data from Short Baseline Experiments 

A combined analysis of Bugey, CHOOZ, CCFR84, CDHS, KARMEN, NOMAD, and LSND for 
the case of a 3+2 mass hierarchy gives two possible solutions ^0]. These are shown in Tabled 
Column 2 of Tabled shows the global best-fit, while column 3 gives the solution if all neutrino 
masses were restricted to lie in the sub-eV range ^. Note that the values of 6*34, ^35 and ^45 are 
not constrained by any experiments to date. However in what follows, we shall assume that all 
concerned mixing angles (including ^13) are large enough so that all the transition probabilities 
are adiabatic for the static density profile of the supernova in absence of shock effects. We shall 
briefiy discuss in section 4.4 the effect of reducing the mixing angles and hence the degree of 
adiabaticity on the resultant supernova neutrino signal. 

3 Supernova Neutrinos and Shock waves 

3.1 Time Dependent Luminosity and Energy Spectra. 

About 3 X 10^^ ergs of energy is released in a type-II supernova, 99% of which is in the form of 
neutrinos. These neutrinos drift out from the dense core of the proto-neutron star and beyond 
a certain radius determined by the energy of the neutrino, travel freely and escape. This radius 
of last scattering which characterizes the energy distribution of the neutrinos is generally called 
the "neutrinosphere" . Neutrinos are mainly released in two phases. A short "neutrinonization 
burst" of pure z/g is produced by electron capture on protons when shock wave crosses the 
neutrinosphere. Subsequently the majority of the neutrinos are produced as u — u pairs of all 
three flavours over a period of 10-20 seconds as the proto-neutron star cools. 

Even if all the six species of neutrinos are dominantly produced by the same mechanism, their 
luminosities and energies at their respective neutrinospheres can differ due the difference in their 
degree of interaction with matter. If the effect of weak magnetism and muon production inside 
the supernova core is neglected, the spectra of z/^, I'r and Ur are approximately equal, and so 
they are often grouped together and collectively denoted as u^. An accurate prediction of the 
time dependent energy spectra of the emitted neutrinos requires a full simulation including all 
significant neutrino interactions. To date only the Lawrence Livermore group [22] has published 
detailed results of the energy spectra of the neutrinos over the full duration of the supernova, in 

^Note that for the best-fit solution shown in column 2, the sum of neutrino masses are in conflict with the 
current cosmological bounds (see for instance ||30| and reference therein). However, we still use them in this 
paper as an illustrative example for the case where the sterile neutrinos are non-degenerate. 
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a simulation with a successful explosion. The resulting luminosity and average energies of the 
neutrinos taken from [22] are shown in Fig. |21 It may be seen that the luminosities of all the 
six types of neutrinos are very similar (apart from the v^. neutronization burst), but the average 
energies of are considerably larger than that of v^. and throughout the duration of the 
supernova. Since neutrino oscillations effectively flip the energy spectra of the and/or v^. with 
that of the and since the extent of this flip depends on the value of the mixing angle and the 
neutrino mass hierarchy, this difference in the initial energy spectra of the different (anti)neutrino 
species has been used in previous studies in the context of three [12] and four jT^I neutrino 
oscillations to put bounds on the neutrino mass spectrum and mixing angles by exploiting the 
energy dependence of the cross sections relevant to the detectors on earth. However, more recent 
work [22] on neutrino transport inside the supernova has cast doubt on the magnitude of this 
effect. When all neutrino scattering processes, some of which were not included in the Livermore 
study, are included, the average energies of the different (anti)neutrino flavors, particularly v^, 
and Vx become very similar. The condition of equipartition of luminosity between the different 
species also breaks down and in some simulations one might flnd the luminosity to be almost 
double that of v^. 

In what follows, we shall use the time dependence of the luminosity as given by the Liv- 
ermore group in our estimates of neutrino oscillation effects as this should not be so sensitive 
to the additional scattering processes not included in the Livermore analysis. However when 
considering the oscillations among the three active neutrinos we will estimate the uncertainties 
by computing the expected signal both for the Livermore luminosities and energy distributions 
and for the case where the energy distributions are equal but the luminosities differ in the man- 
ner just discussed. For the case of oscillation into sterile neutrinos the results are relatively 
insensitive to these differences as the signal does not depend on a difference between and V^. 




Time [sec] 



Figure 3: The time evolution of the luminosity and average energy as calculated by the Lawrence 
Livermore group. 
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3.2 Supernovae Shock Wave(s) 



In this analysis we are particularly interested in the shock waves that form within the supernovae. 
It is believed that when the core of the collapsing star reaches nuclear density, the collapse 
rebounds forming a strong outward shock. This is stalled and then regenerated by a neutrino 
driven wind. During this process both a forward and a reverse shock may form In order 
to get a realistic estimate of the density profile of the shock wave it is necessary to use results 
from numerical simulations. However, since detailed numerical results from such simulations are 
not available to us, we will use a simplified profile for the shockwave used in fT^ l2Uj. In order 
to exhibit the impact of both the forward and reverse shock on the resultant neutrino signal in 
earth-bound detectors, we will consider scenarios where we (i) neglect the effect of the shock 
wave, (ii) consider the effect of the forward shock alone and (iii) consider both the forward and 
reverse shock in calculating the oscillation probabilities. The snap shot of the density profile 
at a certain post-bounce time t for the case of a forward shock only is shown in Fig. and 
the case of both a forward and reverse shock is shown in Fig. Et. The shape of the shockwave 
remains roughly constant in time but the height varies as the shockwave propagates through the 
supernovae (c.f.Fig 1 [18j). The effect of the shock wave on the neutrino oscillation probabilities 
has been discussed before in [T^ ITHj for the forward shock and in 120] for the forward and 
reverse shock. We will briefly review it again in the next subsections for completeness. 
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Figure 4: (a) The density profile of the forward shock as a function of the distance from the 
core of the supernova, (b) The resulting flip probability as a function of time. 



3.3 Effect of the Shock Wave 
3.3.1 Forward Shock Only 

We will consider here the effect of the shockwave on the neutrinos emerging from the supernova 
in the case where all the relevant mixing angles are sufficiently large so that the resonance for 
the static density profile in the absence of the shockwave is completely adiabatic. 

The density structure shown in Fig. 0^ exhibits three separate regions. Initially, for times 
t < tl, the resonant density pji (= (Am^ cos 26')/(2^/2Gi7A^^Y'e£')) lies below the density at 
which the shock wave appears, pa < pi, i.e. in region A. The resonance crossing is adiabatic 
to a good approximation in this region and therefore the flip probability vanishes for t < ti 
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Figure 5: (a) The density profile of tlie forward and reverse shock as a function of the distance 
from the core of the supernova, (b) The resulting flip probability as a function of time. 

as shown in Fig. HJd. As the shockwave moves out pi decreases and for times t\ < t < the 
resonant density lies in region B , pi < < p4^. Due to the large density gradient at the shock 
front, in this region the flip probability increases to cos^ ^ ^ 1 (c.f. eq®) and therefore the 
resonance becomes non-adiabatic in the period ti < t < t4 as shown in Fig. As may be 
seen from the figure, the neutrino pass through two adiabatic resonances and one non-adiabatic 
resonance. The flip probability passing through each adiabatic resonances is approximately 
zero and therefore the total flip probability is equal to the flip probability passing through 
the non-adiabatic resonance alone. Finally in region C the shockwave has passed through the 
resonant density so that p4 < pr and the resonance is adiabatic, therefore the flip probability is 
approximately zero again for t > t^. 

3.3.2 Forward and Reverse Shock 

The situation is illustrated in Fig|3| For the density in the appropriate range there will now be 
two non-adiabatic resonances, corresponding to the points ri and r2. As before at early times 
t < tl, the resonant density pn lies below the density at which the shock wave appears, pr < pi, 
i.e. in region A. As before, to a good approximation in this region the resonance crossing is 
adiabatic and so the flip probability vanishes for t < ti as is shown in Fig. EJd. As the shockwave 
moves out pi decreases and for times ti < t <t2 the resonant density hes in region Bi such that 
pi < Pr < p2- Thus in region Bi only the forward shock has reached the resonant density and 
is therefore equivalent to region B in the forward shock only case. At later times t2 < t < 
the resonant density lies in region B2 where P2 < Pr < P3- In this region the neutrino passes 
through the resonant density at both the forward and reverse shock and both are non-adiabatic 
resonances. Therefore the mass eigenstates are flipped at the reverse shock and then flipped back 
at the forward shock. When the mixing angle is non-zero the flip probability at the forward and 
reverse shocks is cos^ 9 and the total flip probability is | sin^ 29. Between < t < t^, the shock 
profile corresponds to region B^ where ps < Pr < P4 and the neutrinos cross their resonance 
density through the forward shock only. Hence the flip probability again increases suddenly to 
cos^ 9. For t > t4^, we are in region C where p4 < pr and the shockwave has passed through the 
resonant density. The resonance is therefore again adiabatic. The total flip probability of the 
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Model 


{ED 


{ED 


{ED 






LL 


12 


15 


24 


2.0 


1.6 


Gl 


12 


15 


18 


0.8 


0.8 


G2 


12 


15 


15 


0.5 


0.5 



Table 2: The average energies and total fluxes characterising the primary neutrino spectra 
produced inside the supernova. The numbers obtained in the Lawrence Livermore simulations 
are denoted as LL, while those obtained by the Garching group are denoted as Gl and G2. 

neutrinos for the case where we have both forward and reverse shock has the form of a double 
peak as is shown in Fig. El The same structure of the flip probability can be caused by a density 
proflle where the two non-adiabatic resonances which are not "overlapping", such as Fig. IHl 
where the flip probability is only non-zero in regions Di and D^. However this also requires that 
there exists a forward and reverse shock. 
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Figure 6: An alternate density profile giving the same flip probability as shown in flgureEl 



4 Probing Supernovae Neutrino Oscillations in IceCube 

In this Section we present the neutrino induced signal expected in the IceCube detector for a 
galactic supernova event. We first give the results for the case of three active neutrinos and 
then take up the case where two additional sterile neutrinos are added to the neutrino mass 
spectrum. We first turn to the discussion of the properties of the IceCube detector. 

4.1 IceCube as a Supernova Detector 

IceCube is a future km^ neutrino telescope, under construction in Antartica [T7]. On completion, 
it will contain 4800 optical modules deployed into the ice which will detect the Cerenkov photons. 
Though designed to observe ultra high energy neutrinos, IceCube can detect a supernova through 
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its neutrinos by detecting Cerenkov light from e"*" produced by z/g capture on protons. Each event 
cannot be distinguished in the detector, but the addition to the background photon halo coming 
from the Cerenkov radiation associated with the produced positron can be measured. The 
number of additional photons is 

= M^l^^^'^^-^^^^^^^i^ / dEa{E)EF{E,t) (23) 

where fch, fabs and foM are fudge factors of order 1 and are given in the appendix, p = 6.18 x 10^^ 
m~^ is the density of targets in ice, a (in cm^) is the cross section for z/gP — > ne+, E (in MeV) 
is energy, D (in cm^) is the distance to the supernova, t (in s) is time, F (in MeV"'^ s~^) is the 
flux of anti-electron neutrinos. The initial spectra of neutrino species from a supernova is 
parameterized as [HI 

n, , $(t) (a(t) + l)°W+V E Y^'^ f , . E \ 

where {E) and $ are the average energy and total number flux and a is a dimensionless parameter 
which typically takes the values 2.5-5. For the results presented in this paper, we have assumed 
ttg = 3 and = 4. In order to compare the impact of the uncertainties on the average energies 
and fluxes of the neutrinos obtained in different supernova computer simulations, we will present 
our results using supernova neutrino parameters given by both the Lawrence Livermore and 
Garching groups. Specifically, we consider three cases shown in Table |21 . 
The neutrino flux in the detector is 

Fp = Y,F^Pap (25) 

a 

where 

i 

Pa^=J2\Ua,\'P^^ (27) 

j 

P., = \{u.WT)\' (28) 

is the probability that a Iz/^) produced inside the supernova emerges in the i^^ mass 
eigenstate, |z/j). f/™- is the mixing matrix at the point of production, the neutrinosphere. Pij 
is the probability that a {uj^) mass eigenstate in matter appears as |z/j) mass eigenstate in a 
vacuum, this is known as the flip probability. P-® is the probability of detecting the {ui) mass 
eigenstate in the {up) weak interaction eigenstate. In this paper, we do not take into account 
Earth-matter effects and therefore 

P^®P = \U^(^\' (29) 
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4.2 IceCube Background Signal 



IceCube has an irreducible background photon halo with energies comparable to those produced 
by the supernova neutrinos. In the analysis of a supernova signal the additional number of 
photons need to be distinguished from this background. The mean number of detected photons 
is 

rii = NuVlpehin + Ndet (30) 

where Nm is the number of optical modules, vipe is the mean background detection rate per 
optical module, tun is length of the time bin and Ndet is the signal. Using Poisson statistics the 
fluctuations in the background photon halo is given by s/fU.. The background rate is SOOifz per 
optical module and IceCube would have 4800 optical modules. Thus for time bin of 100ms the 
random fluctuations are expected to be a = \/2A x 10^ + Ndet- 

In addition there is an error due to the uncertainty in the time resolution, of the detector 
giving 

An+ = -^(Uj+i +ni_i) (31) 

Ifem 

Anr = —2ni. (32) 

thin 

where thin is the size of each time bin and we take a detector time resolution of 5t = 10 nsec for 
our estimate of the uncertainties. In our analysis we will use as an upper and lower estimate 
of the true event number the values rii + Anf + a and rii — An^ — a respectively, where rii 
is the number detected in the i^'^ bin. In our numerical plots we will show only the additional 
photons produced by the supernova neutrinos, but for the fluctuations we will include all errors 
as discussed above. 



4.3 Results 

4.3.1 Three Active Neutrinos 

We first consider the case of a "standard" supernova at 10 kpc from earth with average energy 
and luminosity for the neutrinos taken from the Lawrence Livermore group simulation. The 
initial spectrum of neutrinos used is that given in Eq. ()24|) . At the neutrinosphere the matter- 
induced potential is much greater than the vacuum mass squared differences; for anti-neutrinos 
this is large and negative. Therefore Vf. is approximately equal to the mass eigenstate with the 
lowest mass squared. In the normal hierarchy this is vi and in the inverted hierarchy this is z/3. 
Fig. ^shows the evolution with density of the mass squared of the anti-neutrino mass eigenstates 
for the normal and inverted hierarchies. From this one may readily see if there is a level crossing 
and determine the flux in the detector using Eqs. (j25|) and (j26p . The relevant expressions for 
the probabilities are given in Table IHl The resonant oscillation effects for anti-neutrinos occur 
only if there is an inverted mass hierarchy. In this case the flux of electron anti neutrinos in the 
detector is given by 

F, = F0 + P,,(F°-F°) (33) 

where 

Pee = PlzWel? + (1 " PlzWeZ? (34) 
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Table 3: The probabilities ^^"^ f*^^ three active neutrinos, where P™' is given by Eq. 
fj26|) and Pij is the flip probability at the resonance between the z/j and Vj mass eigenstates. Only 
Pi3 appears in the expression for the probabilities since P12 = for the current values of Am2i 
and sin^ 6*12. 

From Eq. ()33|) we see that if the initial flux F^{E) of the and Ux are identical, there will 
be no observable oscillation effect. In the LL case the average energy of the I'x is much larger 
than that of the Ug. Since the detection cross-section increases quadratically with anti- neutrino 
energy, flavor oscillations inside the supernova are expected to enhance the anti-neutrino signal 
in IceCube. In Fig. [7|the number of photons that would be detected by IceCube due to the 
supernova anti-neutrinos is shown for the normal hierarchy and for the inverted hierarchy for 
the case of no shock wave, a forward shock wave and a forward and reverse shock. In the figure 
the width of the line represents the upper and lower event number as discussed in the previous 
Section. Thus signals much larger than this width should be observable. 

Due to the energy dependence of the interaction cross section in the detector the more 
energetic component will give the larger signal. In the case of a normal hierarchy there are no 
resonant level crossings for the antineutrinos and the flux on earth originated as the Ue flux 
at the neutrinosphere (c.f. Table El). Since in this case the resultant neutrino flux has the least 
energetic component, it explains why the normal hierarchy case in Fig. [7| corresponds to the 
lowest bound for the number of photons in the detector at all times. For inverted hierarchy 
if we ignore the effect of shock waves, then for large values of ^13 the transition is completely 
adiabatic and the resultant Ue flux on earth originated as the Ux flux at the neutrinosphere. Since 
in this case the z/g flux at the detector is made up entirely of the most energetic component, 
the resulting number of photons provide the upper bound on the expected number of photons 
in the detector. The effect of the shock wave is to bring about abrupt changes in the oscillation 
probability as discussed in the Section 13.3.21 In particular, we had seen that the transition 
probability flips from being completely adiabatic to almost non- adiabatic at the shock front. 
Therefore, including the effect of the shock wave means that the Ue flux on earth at the time 
when the shock front crosses the resonance region originated as the flux at the neutrinosphere, 
causing the reduction seen in the signal at these times in Fig. [7| We see that the structure in 
the signal for the forward and reverse shock does show a single and double bump structure 
respectively, conforming to the shape of the oscillation probability discussed in Section 13.3.21 

In both the single and double shock cases the structure in the signal above the photon 
background is clearly visible, demonstrating that IceCube will indeed be capable of distinguishing 
between the normal and inverted hierarchy for the case the initial energy spectra of different 
flavours differ. However, as discussed in Section IXTl the inclusion of further scattering processes 
within the supernovae is expected to have the effect of reducing the difference in the average 



16 



□ 


13 


N3 




13 


FRS 




13 


FS 




N3 







4 6 S 10 

Time (s) 



12 



Figure 7: The left hand panel shows the number of photons (A^) that would be detected in 
IceCube for the case of three active neutrinos only with either normal or inverted hierarchy. 
For inverted hierarchy we show the results where we disregard the shock effect (NS), consider 
only the effect of the forward shock (FS) and take effect of both the forward and reverse shocks 
(FRS). For normal hierarchy, the shock does not have any effect and we have the results denoted 
as N3 in the figure for all cases. The right hand panel shows the difference (AA^) between the 
number of photons expected in presence of shock to the number of photons in absence of shock 
effects. The observed photons are split into 100ms time bins and the width of the lines reflect 
the expected fluctuation in the photons detected. 



energy between the z/g and Vx components. In this case the structures in the signal just discussed 
and indeed the effect of oscillations itself due to the energy differences between the different 
(anti)neutrino flavors will not be present. However all is not lost because the effect of the 
additional scattering processes, while reducing the energy difference, increases the difference in 
number flux of the neutrinos with the number flux of the Ve becoming almost twice that of the 
Vx components during the accretion phase (cf. Table E))- To make this quantitative we have 
calculated the number of photons expected in the IceCube detector for the cases Gl and G2 of 
Table |21 Since results for the Garching simulations are not available over the full time period 
of the supernova, we have used the average energy and total fluxes for Gl and G2 from Table 
121 but the time dependent fall of the flux from the simulations of LL. The results are presented 
in Fig. IHl where the upper panels are for the case Gl and lower ones for G2. The results we 
get are very similar to the one we had for the LL case. The reason is that even though the 
ratio of {E^J/ {E^^) is nearly 1 for the Garching simulations, the ratio of the fluxes $p^/<I>|^^ is 
nearly half. This means that in this case even though neutrino oscillation does not change the 
energy of the resultant z/g spectra emerging from the supernova, it changes the total number 
flux, increasing the number of Ve and hence producing an enhanced signal in the detector. 

Note that the mere existence of structure due to the shock waves provides the evidence for 
the inverted hierarchy. However as is evident from Figs. [7| and |H1 the signal contains much 
more information, the sign of the effect giving information about the relative importance of 
an energy difference or luminosity difference between the z/g and components. Also the 
shape of the signal provides seismological information capable of determining much about the 
nature of the supernova shock wave or waves. Clearly, given the sensitivity of the signal to the 
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Figure 8: Same as Fig. [3 but with equal initial energy spectra for the active antineutrinos and 
different luminosities as discussed in the text, (a) and (b) are for the case called Gl while (c) 
and (d) show the case G2. 

initial spectra and luminosities, it is desirable to have a full time dependent simulation of the 
luminosity and energy spectrum of the Lawrence Livermore type but including all significant 
neutrino interactions within the neutrinosphere. 

4.3.2 Three Active and Two Sterile Neutrinos 

We turn now to the case in which there are additional sterile neutrinos capable of explaining 
the LSND anomaly. Once again resonant effects due to the passage of shock waves through 
the supernova can give significant information about the neutrino mixing. In this case we do 
not have to rely on a detailed knowledge of the luminosities and energy spectra of the active 
neutrinos because oscillation to a sterile neutrino necessarily corresponds to a reduction in the 
signal and therefore it's much easier to observe oscillation and shock effects in this case. 

The analysis proceeds in a manner similar to that for the three neutrino case. We start with 
the evolution of the mass squared of the mass eigenstates for the different possible hierarchies 
as shown in Fig. |21 From this it is straightforward to compute the v^. flux in the detector using 
Eqs. f03|) and Table El gives the probabilities P"^ and for the different neutrino mass 

spectra. This simplifies for the case considered here, |?7eiP, |f^e2p >> |f^e3p, |^e4p, |^e5p giving 
the results in Table|S| The angles ^14 and ^15 are taken from the best fit values of tabled We take 
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Table 4: The probabihties P^ and P™ for three active neutrinos plus two sterile neutrinos, 
where P^ is given by Eq. and Pjj is the flip probability at the resonance between the Ui 
and Uj mass eigenstates. 
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Table 5: The flux of neutrinos in the approximation that jf/eiP, \ Ue2\ 
and Pii, ~ P2i, where i=4 or 5. 
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634 = 635 = 10~^; note that these are not constrained from any data. Finally we take 613 = 0.142, 
within the bounds from equation l2.3.3l and 624 = 0.12 and 625 = 10~^, chosen such that the LSND 
mixing angle Olsnd is equal to it's best fit value, where sin^(26'L5Ar£)) = 4:{Ue4U^4 + Ue^Ufj,^)'^. 
These angles are large such that the transition probabilities change from being adiabatic to 
non-adiabatic as the shock wave passes through each resonance. The signal is dominated by 
detection of anti-neutrinos leaving the supernova in the first and second mass eigenstate. In 
all but one case the mixing angles are sufficiently small such that the multiple resonances that 
these anti-neutrinos cross are far enough apart such that they do not overlap and therefore each 
resonance can be considered as independent. The exception is the combination of resonances 
corresponding to Am^g and Aml2, however numerical checks have shown that the total transition 
probability is still well represented by a product of the two flip probabilities. Using this we have 
computed the number of photons that would be detected in IceCube for a galactic supernova. 
Figs, ini GDI HH UH UH and [T31 are for the LL case while Fig. [T31 compares the impact of the 
change in the initial neutrino fluxes and spectra from LL to Gl and G2. 
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Figure 9: The number of photons detected for a N2+I3 mass hierarchy calculated for 100ms 
time bins with a forward shock (FS) and without a shock (NS). 
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Fig. El shows the number of photons detected in time bins of 100 ms, due to a supernova for 
neutrinos in a N2+I3 mass hierarchy, with and without a shock wave. There are bumps in the 
signal due to the propagation of the shock wave in the first 2s and from 4s-lls. As has been 
discussed for the case of just three active neutrinos, the latter structure is due to the propagation 
of the Shockwave through the resonant density corresponding to the atmospheric mass splitting. 
The bumps in the first 1.5s is due to the propagation of the shockwave through the resonant 
densities corresponding to the sterile neutrinos. 

We first concentrate on the Amg^ driven "atmospheric resonance" to determine the changes 
due to the mixing with sterile neutrinos. The effect of the shock wave is only relevant if the 
three active neutrinos have an inverse hierarchy in their masses. This is evident from Table 
where the flip probability P13 appears only for cases with 13. For these cases we show in Fig. 
ITUl the number of photons detected on the same time bins as that shown in Fig. [3 We do 
not show the case for H2+I3 because, c.f. Table El the shock wave only causes a change in the 
signal if it also passes through a sterile resonance at the same time. In theoretical simulations 
of the shockwave this does not happen and therefore no "bump" should be observed. From the 
Figs. [7|and^lwe see that the N2+I3 hierarchy gives an increase in the signal compared to the 
case of only 3 active neutrinos while the 12+13 case shows a decrease. In both cases we get a 
statistically significant signal. 

Turning to the sterile resonance region we present in Figs. ^2 El El and our numerical 
results for the sterile resonances. The structure of the signal spans of the order of 0.1s and 
therefore a detector needs to have a time resolution of at least 10ms to try to disentangle the 
information. The number of bumps is determined by several factors. First, it is dependent 
on the mass hierarchy which determines the number of resonances which the neutrino passes 
through, (see Table E}. If a neutrino passes through two sterile resonances it could create two 
bumps; however if the resonances are close inside the supernova, only one broader bump may 
be observed. If the shockwave has a forward and reverse shock the double peak structure in the 
probability can cause several bumps in the signal, as can clearly be seen in Fig. El(c) and (d). 
For the N2+I3 and H2+N3 mass hierarchies, because |f/eip, |t42p >> l^esP, |f^e4p) l^^sP, the 
number of events in the absence of a shock is an order of magnitude smaller than that of other 
hierarchies. Therefore observation of a large ratio of the height of the bump to the background 
supernova signal would be characteristic of these mass hierarchies. However in practice this 
could be misidentified if the luminosity varied rapidly in time. 

The effect of changing the average energy and total number flux of the neutrinos is shown in 
Fig. El The six panels are for the six different mass spectra and we consider the case where we 
have both forward and reverse shocks. We can see from the figure that changing the supernova 
neutrino model from LL to Gl or G2 does not wipe out any of the interesting structures in the 
signal. 

5 Summary and Conclusions 

Our analysis of the neutrino signals in the IceCube detector that would result from a "standard" 
supernova at 10 kpc from earth has shown that oscillations to sterile neutrinos will be easily 
detectable for a wide variety of neutrino spectra. As may be seen from the figures, for 10 msec 
time bins, the event rate is typically very large and supernovae as far as about 30 kpc from 
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Figure 10: The number of additional photons detected in 100ms time bins, for the cases of 
a forward shock (FS), a forward and reverse shock (FRS) and no shock (NS) for the mass 
hierarchies: (a) N2+I3, (b) 12+13. The lower panels show the corresponding difference in between 
the number of photons expected due to the effect of the shock for the mass hierarchies: (c) N2+I3, 
(d) 12+13. 

earth will still give observable effects. The reason the signals are so distinctive is a result of 
the combination of resonant conversion within the supernova with the time dependent structure 
resulting from the propagation of a shock wave through the supernova. Together this gives a 
characteristic time dependent signal that is readily observable with a detector time resolution 
better than 10 msec. Since oscillation to sterile neutrinos generates a disappearance event the 
results are relatively insensitive to uncertainties in the initial energy spectra and luminosities. 

The signal coming from resonant effects involving just three antineutrinos with an inverted 
mass hierarchy is much more sensitive to the initial spectra and luminosities because the signal 
vanishes if both are degenerate for the three flavours. For the case the energy spectra differ 
significantly, as in the Lawrence Livermore simulation, the resulting signal is again clearly visible 
in IceCube. However if the energy spectra are degenerate, as recent studies seem to indicate, 
then one must rely on differences in the initial himinosities. To date no full simulation has 
been performed over the time period of interest but extrapolating the recent simulations, which 
extend up to 750 ms, to later times suggests that even in this case the resonant conversion in the 
presence of a shock wave would be clearly visible despite a reduction of about a factor of two in 
the signal. To improve our confidence in this result it is clearly of importance to have a reliable 
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Figure 11: The number of photons expected in 10ms time bins, due to a forward shock (FS) and 
a forward and reverse shock (FRS), for the mass hierarchies: (a) N2+N3, (b) N2+I3, (c) H2+N3, 
(d) H2+I3, (e) I2+N3. We show the results for both the solutions given in Table [TJ Note that 
the shock wave has no observable effect for 12+13. 
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Figure 12: The difference in tfie number of expected photons in 10ms time bins due to the effect 
of a forward shock (FS) and a forward and reverse shock (FRS), for the mass hierarchies: (a) 
N2+N3, (b) N2+I3, (c) H2+N3, (d) H2+I3, (e) I2+N3. 
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Figure 13: The number of photons detected in 10ms time bins due to the Shockwave in various 

mass hierarchies, where (a) is for forward shock and Am^^ = 22eV'^, Am^^ = 0.92eV'^, (b) is for 
forward shock and Amg^^ = 0.89eV'^, Aml^ = OAGeV'^, (c) is for forward and reverse shock and 
Am|i = 22eV^, Aml^ = 0.92eV^, (d) is for forward and reverse shock and Am|i = 0.89eV^, 
Ami = OAGeV^. 
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Figure 14: The number of additional photons detected in 10ms time bins due to the Shockwave 

in various mass hierarchies, where (a) is for forward shock and Am^^ = 22eV'^, Am^^ = 0.92eV^, 

= 0.89eV'^, Am|-^ = OAGeV^, (c) is for forward and reverse 
22eV^, Am|i = 0.92eV^, (d) is for forward and reverse shock and Am|i ~ 
O.SQeV^, Am^i = 0.466^^. 
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Figure 15: Effect of clianging tlie supernova neutrino model on tfie resultant signal in the 
detector. The number of photons expected in 10ms time bins, due to a forward shock and 
Aml^ = 0.89eV^, Am|i = 0.466^^, for the mass hierarchies: (a) N2+N3, (b) N2+I3, (c) 
H2+N3, (d) H2+I3, (e) I2+N3 (f) 12+13. 
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simulation of neutrino propagation through the supernova over a 10 second period including all 
the significant antineutrino processes. 

Observation of sharp time dependent neutrino oscillation signals will also provide a probe of 
the propagation of the shock wave through the supernova. As our simulations show it is possible 
that one can see both the forward and reverse shocks and the structure of the signal will give 
information about their spatial structure. At present our knowledge about the nature of the 
shock is quite primitive and clearly better simulations of the shock properties will improve our 
ability to interpret the signals. A very recent paper [T^ has pointed out that turbulence after 
the shock can significantly affect the signals, broadening the structure, and clearly it is also of 
importance to study this in more detail. 

In summary we have found that the detailed time structure of the neutrinos coming from 
supernovae can provide information, not accessible by other means, both on the mass spectrum 
and mixing angles of neutrinos and on supernovae seismology. This applies to those neutrino 
mass hierarchies and mixing angles, both for three active neutrinos and for a combination of 
active and sterile neutrinos, which give rise to resonant transitions within the supernovae. For 
the case involving sterile neutrinos our analysis has concentrated on the parameter space which 
can explain all neutrino oscillation phenomena including the LSND measurement. However, 
even if the LSND result should not be confirmed, the sensitivity of experiments such as IceCube 
is such that they will still be able to probe for a sterile neutrino component not otherwise 
observable. 
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Appendix: Detector properties 

The following analysis is the same as that presented by Dighe, Keil and Raffelt . For simplicity 
only events detected by the v^P ~^ f^^^ are considered. Each e"*" emits Cerenkov light according 
to 

d'^N^ _ 2na sin^ 9 ^^^^ 

where A'^ is the number of emitted photons x is the distance traveled by the e"*", A is the 
wavelength of the emitted photons, a is the fine structure constant and 6 is the angle of the 
emitted photons to the direction of the e"*" such that cos 6' = {n(3)~^, where n is the refractive 
index of the ice (njce = 1.31) and P is the ratio of the speed of the positron in the medium to 
that in a vacuum {f3ice ~ 1)- Integrating over the observable wavelengths, taking the mean free 
path of a positron to be 12 cm for an energy of 20 MeV and assuming the two are proportional 
the number of photons emitted between observable wavelengths is: 

^ = 191MeV-'fch (36) 
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where E^^ is the energy of the and fch is a fudge factor. The number of detected photons 
per optical module is: 



_ 47.75 



1 rtmax roo 

fchfabsfoMP^ dt dEa{E)F{E,t)E (37) 

_r Q -^cat ^acc (QQ\ 

^^^-0:20250^1^ ^^^^ 



fobs = i?a6./100m (39) 

where p = 6.18 x 10^^ m~^ is the density of targets in ice, a (in cm^) is the cross section for 
T>eP ne~^, E (in MeV) is energy, D (cm^) is the distance to the supernova, t (in s) is time, F 
(in MeV~^ s~^) is the flux of anti-electron neutrinos, Q is the average quantum efficiency of the 
detector, Acat is the effective photo cathode detection area, flacc is the angular acceptance range 
and Rabs is the absorption length of photons by the ice. 
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